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We discuss the successes of the dynamical mean field (DMF) approach to metal 
insulator transitions in both the clean and the disordered limit. In the latter 
case, standard DMF equations are generalized in order to incorporate both the 
physics of strong correlation and Anderson localization effects. The results give 
new insights into the puzzling features of doped semiconductors. 



1. Introduction 

How a substance evolves from a metallic to a non metallic state is one of the 
most fundamental and richest problems in condensed matter physics. 

In general, there are several mechanisms at play. Electron-electron interactions 
can drive a metal to insulator transition in a pure substance (Mott, 1990). This 
transition is named after Sir Nevill Mott who laid down the foundations for the 
physical understanding of this phenomenon. Another route was discovered by 
Anderson (1958) who realized that sufficiently strong disorder can drive a metal- 
insulator transition even in systems of non interacting electrons. The theoretical 
description of the situation when both effects are present is a central unsolved 
problem. For recent reviews see Lee & Ramakrishnan (1983) and Belitz & Kirk- 
patrick (1995). 

Early treatments (Abrahams et al. 1979, Wegner 1976, Schaffer &: Wegner 
1980, Finkelshtein 1983, 1984) used direct analogies from the theory of magnetism 
and the scaling approach to critical phenomena. The technical apparatus of the 
approach was a field theoretical non linear sigma model and an expansion near 
two dimensions. The physical content of such theories was expressed concisely 
in terms of an extension (Castellani et al. 1987) of the Fermi liquid approach 
to disordered systems. These ideas were very successful in the description of 
the transport and thermodynamic properties of weakly disordered metals. The 
approach, however, encounter several difficulties in accounting for experimental 
observations in many systems. The origin of these difficulties can be traced to (1) 
strong interaction effects, and (2) strong statistical fluctuations. One important 
manifestation of this effect was the formation of local magnetic moments, objects 
which can not be understood from a weak correlation perspective. This aspect 
of the problem and its theoretical treatment were addressed in Paalanen et al. 
(1988) and Milovanovic et al. (1989). 
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In the last few years, a new dynamical mean field approach to the strong 
correlation problem has been developed (Metzner & Vollhardt f989, Georges & 
Kotliar f992). It has explained many puzzling features of clean three dimensional 
transition metal oxides (for a review see Georges et. al. 1996). Very recently, it 
was extended (Dobrosavljevic & Kotliar 1997) to incorporate the interplay of 
Anderson localization and correlations effects. 

The mean field method is largely complementary to the scaling approach. While 
the latter concentrates on the long wavelength modes which presumably govern 
the immediate vicinity of the transition, the former focuses on the local charge 
and spin dynamics of the electrons. The latter is based on an expansion around 
the lower critical dimension, the former being formally valid in large dimensions. 
It is our view that a good mean field understanding of the metal to insulator 
transition problem is a necessary first step towards a comprehensive theory. 

In this talk we will not focus on the technical aspects of the method which are 
thoroughly reviewed elsewhere. Instead we will attempt to give a simple descrip- 
tion of the physical content of the approach, and of the results obtained by this 
method. Section |2| introduces the essential idea of the mean field method for a 
general periodic solid. It is our hope that its generality and local character will be 
of interested to the chemistry community which has stressed the importance of 
local bonding and correlations through the years. Section [| discussed the exten- 
sion to the disordered case. We stress the importance of statistical fluctuations 
and how they are captured by the mean field method. In section || we discuss 
the physics of the density driven Mott transition and explain why the mean field 
approach works so well for three dimensional compounds, stressing the role of or- 
bital degeneracy. We then move to the physics of disordered interacting systems 
which was the focus of many years of Mott's research. Section || discusses our 
views on the puzzle posed by the differences between uncompensated semicon- 
ductors and disordered alloys. In section 5, ^ and 7 we elaborate on the physical 
content of the mean field approach and argue it captures some essential elements 
needed to resolve this puzzle. We conclude in section ||] with a comparison to 
other approaches and outline some directions for further work. 



2. Dynamical mean-field theory (DMF) theory - the clean limit 

To motivate the dynamical mean field (DMF) approach, it is useful to draw 
some analogies with a well established approach to the electronic structure prob- 
lem, the density functional theory (DFT). The basic physical quantity in DFT 
is the density p(r), and the free energy of a system is written as a functional of 
the density -FdftIpC'")]- Its minimum gives the physical density of the system in 
question. In practice, the form of the functional Fdft is unknown and various 
approximations such as the LDA (local density approximation) are used instead 
of the exact density functional. 

The dynamical mean field approach is very similar in spirit, except that it 
adopts the local spectral function A(oj, r) = —^ImG(r,r;u! + iS) as the basic 
quantity. One then writes the free energy as a functional of this quantity, and 
its minimization leads to the dynamical mean field equations which give the one 
particle spectral function for the system of interest. One can motivate the exten- 
sion from the density to the spectral function (which can be thought of loosely 
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as a kind of energy-resolved density ) in the context of strongly correlated elec- 
tron systems, by thinking about their photoemission (and inverse photoemission) 
spectra which demonstrate the existence of bands which have the character of 
atomic configurations (Hubbard bands), in addition to the ordinary quasiparticle 
bands which are analogous to those of the non interacting system. One can incor- 
porate quasiparticle bands and Hubbard bands naturally in a single theoretical 
framework by resolving in energy, the local density, i.e. the local density is the 
integral of the local spectral function over frequency. By resolving in energy the 
localized and itinerant component of the electron, DMF approach treats coherent 
(quasiparticle like) and incoherent (Hubbard bands like) excitations on the same 
footing. It is a unified framework for the description of localized and itinerant 
electrons. 

To illustrate the generality of the method let us start from a Hamiltonian 
containing several orbitals per unit cell. We use a compact notation where the 
index a = (m, a) combines the orbital m and the spin a. 

^lattice ^ . CjQ, tia, 

(ij) i i 

(2.1) 

We now focus on a single unit cell, and integrate out all degrees of freedom 
except for those which reside in the selected unit cell. These are described by 
operators c a and no longer carry a site index. The dynamics of the resulting 
problem is described by an impurity model which describes an impurity (c a ) 
coupled to a bath of fermions (ctfeu) 

H imp = ^2 ( Ea P ~ V C t C P + X! r a/?7<5 c t c c ^ C S + ^2 e b^ A^i/i 

a/3 i bfi 

+ ]T(v^ a a+c Q + /i.c.) (2.2) 

k 

From the impurity model we can obtain all the local correlation functions, since 
by construction the local lattice Green's functions are identical to the impurity 
Green's function G = (G Qj ^) . We use a matrix notation so that the local Green's 
function is given by 

G a ,p(r - r') = — (T T c a {r)c + p {T')) (2.3) 

This should be viewed as a functional of the parameters and V^q,. To 
determine these parameters we construct the "Weiss field" which describes the 
effect of the rest of the electrons on the selected cell, Gq 1 (iuj n ) = (iu n + /J,) I — 

E - (j2b/j. V lZ-ebf )> and the self ener §y of the im P urit y> S(iw„) = Go\iu; n ) - 
G~ l (iuj n ), viewed as a functional of e;, M and V^. These parameters are determined 
by requiring that the bath and the local degrees of freedom describe the electrons 
in the original lattice problem. Namely we can construct the local Green's function 
from the lattice Green's function obtained by adding a k independent self energy 



to the non interacting lattice Green's function (obtained from 2.1 by setting the 
interaction terms to zero) or from the impurity model. 
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G(iuj n ) = Y<((^n + v)i-i(k)-£(iLu n )-E) 1 (2.4) 
k 

The reader will notice many common features between DMFT and the Bragg 
Williams theory of magnetism. However in fermionic problems a new feature 
emerges-the cavity field acquires nontrivial time dependence, allowing a nonper- 
turbative treatment of local dynamics, which proves to be of crucial importance 
for strongly correlated electrons. In particular, the approach incorporates inco- 
herent (inelastic) processes even on this mean-field level, as opposed to most other 
treatments. As a result, the formulation can be used even in the study of non- 
Fermi liquid metallic phases, for example in extended Hubbard models (Si and 
Kotliar 1993). 

3. Dynamical Mean Field Theory of the Density Driven Mott 
Transition in Pure Systems 

The Mott transition in transition metal oxides has received renewed theoretical 
and experimental attention. On the experimental side, new compounds have been 
synthesized (Tokura et al. 1993, Okimoto et al. 1995) and known compounds such 
as V2O3 and NiSe\_ x S x have been studied. On the theoretical side, new insights 
have been obtained from studying the one band Hubbard model in the limit of 
large lattice coordination. In this limit the pressure and temperature driven Mott 
transition as well as the density driven transition can be thoroughly analyzed. 

Several quantitative comparisons between the physics of three dimensional 
transition metal oxides and the one band Hubbard model have already been 
performed (Kotliar & Kajueter 1996, Rozenberg et al. 1995). For example, the 
doping dependence of the electronic specific heat, the resistivity and the Hall 
coefficient in La x Sri_ x TiOs can be explained by the one band Hubbard model 
without adjustable parameters after the values of U and D have been extracted 
from photoemission data (Kajueter et. al. 1996). Similarly, the single band Hub- 
bard model can describe the temperature dependence of both, the optical and 
the DC-conductivity in V2O3 (Rozenberg et al. 1995). 

The Mott transition with integer occupation has been the subject of extensive 
reviews (Georges et al. 1996). Here we focus on some qualitative aspects of the 
physics of the density driven Mott transition, the emphasis will be on the question: 
why does the mean field theory work so well for three dimensional transition metal 
oxides? 

We start with a brief discussion of the qualitative content of the mean field 
theory of the doped Mott insulator. At low temperatures the mean field theory 
describes a Fermi liquid with a an energy scale (renormalized Fermi energy e^*) 
which is proportional to the doping (i.e. deviation from half filling) 5. In infinite 
dimensions, the renormalized Fermi energy is proportional to the quasiparticle 
residue Z which also vanishes linearly in doping, i.e. 6f* = ZD . In this regime, 
the linear term of the specific heat is inversely proportional to 5 and the Hall 
coefficient is unrenormalized from its band structure value. The Fermi liquid 

3 

description is valid up to a temperature scale T that surprisingly scales as 62, 
and is therefore quite small very close to the Mott transition (Kajueter et. al. 
1996). 
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The qualitative success of the mean field approach when applied to three dimen- 
sional transition metal oxides is due to the orbital degeneracy of three dimensional 
transition metal oxides. Orbital degeneracy causes a close competition between 
ferromagnetic and antiferromagnetic tendencies. The net spin-spin interaction is 
a sum of antiferromagnetic and ferromagnetic terms which tend to largely cancel. 
An illuminating example that can be studied exactly is a two site system was de- 
scribed by Kajueter and Kotliar (1997) mimicking a Ti ion in LaTiO^. The basic 
scale of the problem is a bandwidth which is of order of 1 eV, but the splitting 
between different spin configurations turns out to be a much smaller ( between 
10^3 and 1(T 4 eV). 

If we insist on using a one band model which ignores the orbital degeneracy, to 
reproduce the gross features of transition metal oxides, one must introduce in the 
model terms that suppress the tendency towards magnetic order. In infinite di- 
mension, this can be accomplished by choosing lattices with longer range hopping 
matrix elements which induce magnetic frustration or by adding and additional 
ferromagnetic interactions to the Hamiltonian to compensate for the antiferro- 
magnetic exchange which is inherent to the Hubbard model on the hyper-cubic 
lattice. A Hamiltonian of the form: 



with an additional ferromagnetic interaction J increasing as a function of x, can 
serve as a crude caricature of the La\- x Y x TiO^ system. This compound is fer- 
romagnetic for x near 1 but becomes antiferromangetic at small values of x. 
(Okimoto et. al. 1995). 

While this approach is certainly too crude to describe details of the magnetic 
ordering, it expresses the fact that the localization tendencies in this system are 
not caused by, and are largely independent of, the magnetic ordering (which can 
be either ferromagnetic or antiferromagnetic). In the limit of large dimensions 
the J term does not affect the one particle properties, so the calculations of the 
effective mass perf orm ed in the absence of this term apply to this model as well. 
The Hamiltonian (|3.l| ) was used by A. Georges & L. Laloux (1996) to model 
the Wilson ratio of liquid He^, and by Obermeier et. al. (1996) to model the 
magnetic response of cuprate superconductors. 

To summarize, the mean field theory works very well in systems where the 
charge and spin fluctuations are reasonably local. This is the situation in three 
dimensional transition metal oxides (as a result of orbital degeneracy) and in the 
disordered systems described in the next section. 

While the qualitative success of dynamical mean field theory is easy to un- 
derstand, the surprising quantitative agreement of this approach as applied to 
orbital degenerate systems was clarified only recently( Kotliar &: Kajueter 1996, 
Kajueter & Kotliar 1997). It was shown that: the quantitative value of physical 
quantities such as the effective mass, and the optical gap, near the Mott insu- 
lating state with one electron per site which had previously been compared with 
experiments depend weakly on the band degeneracy (the differences are of less 
than 10 percent). On the other hand the high energy behavior of the spectral 
functions have a sizable dependence on the number of orbitals per site. 
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4. The metal to insulator transition in disordered systems: some 

outstanding puzzles 

The presence of disorder adds a new dimension to the metal to insulator tran- 
sition problem. The continuous nature of the metal to insulator transition was 
first demonstrated in this class of systems ( for a general review of metal insula- 
tor transitions see Mott (1990)). As a result of a series of intensive experimental 
studies a basic picture of the behavior of physical quantities near the transition 
has emerged. In this section we describe briefly some aspects of the experimental 
picture which we regard as well established, but which still remain a challenge 
for the microscopic theory. For more comprehensive reviews of the experimen- 
tal situation in this field see Paalanen & Bhatt (1991) Sarachik (1995) and v. 
Lohneyesen (1997). 

It is well established the T = conductivity vanishes as the transition (i.e. 
critical dopant concentration in doped semiconductors) is approached from the 
metallic side. We leave the immediate vicinity of the metal to insulator transition 
out of this analysis because the extrapolation to zero temperature is problem- 
atic (for a good discussion of this point see v. Lbhneysen (1997)). In the region 
where this extrapolation is unambiguous, the conductivity vanishes in a power 
law fashion, a ~ (n — n c )^. In uncompensated materials (one electron or hole 
per dopant ion), the conductivity exponent takes an anomalously small value 
H « 1/2 (Paalanen et al. 1980). Such a small value of the conductivity expo- 
nent proved to be notoriously difficult to explain by any known theory. On the 
other hand, well compensated doped semiconductors (only partial filling of the 
impurity band) behave very differently; their conductivity vanishes with an expo- 
nent jJL ~ 1. Amorphous alloys behave the same way. In the presence of a strong 
magnetic field the conductivity vanishes with an exponent [i ~ 1, in both com- 
pensated an uncompensated semiconductors (Dai et al. 1991, 1992). The presence 
of spin-orbit (SO) coupling (Dai et al. 1991, 1992) does not seem to be a deter- 
mining factor for the behavior of the conductivity. Both Si:B (where SO coupling 
is strong) and Si:P (where SO coupling is weak) behave qualitatively in the same 
way. Uncompensated doped semiconductors also display an anomalous temper- 
ature dependence. In the metallic phase for concentrations not too close to the 
transition the conductivity increases as the temperature is lowered. This behav- 
ior can be reversed by applying strong magnetic fields, resulting in a decreasing 
conductivity at low T. 

In our view these anomalies occur over a wide range of concentrations and 
are not restricted to a tiny critical region in the vicinity of the transition. For 
example, the anomalously small value of the conductivity exponent fi ~ 1/2 
describes the data all the way to n ~ 4n c ! It is then natural to conclude that the 
observed behavior should not be identified with an asymptotic critical behavior 
associated with a narrow critical region. Rather, it should be described by an 
appropriate equation of state, that is expected to follow from a relevant mean-field 
description of the problem. We also emphasize that the anomalies are associated 
with uncompensated systems, where we expect the effects of electronic correlation 
to be the strongest. It is thus natural to approach the problem from a DMF 
perspective. In this picture, many of the above features, but in particular the 
compensation (filling) dependence follow naturally. 

In contrast to the transport properties, the thermodynamic quantities x an d 
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7 vary smoothly (Paalanen et al. 1988) as a function of concentration across the 
transition, and seem to diverge at low temperatures both on the insulating and 
on the metallic side of the transition. The NMR experiments (Paalanen et al. 
1985, Alloul & Dellouve 1987) portray a strongly inhomogeneous picture. There 
is a wide distributions of Knight shifts on the Phosphorus sites. As the transition 
is approached a large number of sites acquire knight shifts that are larger than 
the measurable range indicating the formation of local moments. The knight shift 
on Si is a smoother function of concentration ( Stoddart et al. 1992) suggesting 
that the metal insulator transition takes place in the phosphorus impurity band. 

Underlying this experimental picture is a broad distribution of energy scales 
which makes the metal-insulator transition in disordered systems very different 
than in the clean case. In addition, this behavior reflects unusually developed spa- 
tial fluctuations of the electronic system. While such behavior does not normally 
occur in ordinary metals, it is a natural consequence of the incipient localization 
of the electronic states - which by definition cannot happen in a uniform fashion. 

5. Statistical DMF theory 

We have argued that in presence of disorder, the situation is qualitatively more 
complex than in the clean limit. Clearly, it is necessary to extend the DMF ideas 
in order to capture the crucial effects of disorder - the spatial fluctuations of the 
order parameter. In the following we present a simple pedagogical derivation of 
the statistical DMF theory, and discuss some of its main features. 

We follow an approach very similar to the Thouless- Anderson-Palmer (1977) 
formulation of the mean field theory of spin glasses. Specifically, we treat the 
correlation aspects of the problem in a dynamical mean-field theory fashion, but 
allow spatial variations of the order parameter in order to allow for Anderson 
localization effects. The theory is then exact in the noninteracting limit, and 
reduces to the standard dynamical mean field theory in absence of disorder. 
The approach can be applied to any lattice model of interacting fermions. For 
simplicity, we consider a simple single-band Hubbard model with random site 
energies given by the Hamiltonian 

ij <y i 

Following the general spirit of the DMF theory, we focus on a particular site 
i of the lattice, and integrate all the other sites. This procedure is formally ex- 
act, but the resulting local effective action takes an arbitrarily complicated form 
(Dobrosavljevic & Kotliar 1994), containing n-point vertices of all orders. How- 
ever, within the DMF approach, one retains only the contributions quadratic in 
Fermion fields, and the local effective action assumes the form 

Seffii) = £ f dT f ^c{ a {r){b{T - tO (d T + e z -v) 

a J o Jo 

+A i)ff (r, t0)cv(t0 + U f drn h1 (T)n id (T). (5.2) 

Jo 

Here, we have used functional integration over Grassmann fields Cj )(r (r) that rep- 
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resent electrons of spin a on site i, and nj )(T (r) = c\ a (T)ci^{T). The "hybridization 
function" Aj(r, r') is given by 

A i (u n )= £ 4-G«K). (5.3) 

j,fc=l 

The sums over j and runs over the z neighbors of the site i, and 

Gf k K) =< c] K)c fc (w n ) >W (5.4) 

are the lattice Green's functions evaluated with the site i removed. In general, 
these objects can be expressed through ordinary lattice Green's function as 

Gf k = G ]k - (5.5) 

We emphasize that the above construction is carried out for a fixed realization of 
disorder defined by a given set of random site energies {ei}. The DMF truncation, 
which keeps only the quadratic contributions to the effective action is exact for: 

(A) infinite coordination (z — ► oo) (Dobrosavljevic & Kotliar 1993, 1994), or 

(B) non- interacting electrons for arbitrary coordination (Anderson 1958). 

Note that in the case (A) the sum in the definition of Aj(u; n ) runs over in- 
finitely many neighboring sites, so that the hybridization function is replaced by 
its average value. As a result, the all the spatial fluctuations in the "cavity" rep- 
resenting the environment of a given site are suppressed, prohibiting Anderson 
localization effects. For noninteracting electrons, the resulting theory reduces to 
the well known "coherent potential approximation" (Elliott et al. 1974). 

Since one of the main goals of the statistical DMF theory is to incorporate 
Anderson localization effects, we concentrate on finite coordination lattices. In 
this case, Aj can be considered to be a functional of the lattice Green's functions 
Gjk, evaluated for fixed disorder, i. e. 

At = A, [G jk ]. (5.6) 
At strong disorder, Aj will exhibit pronounced fluctuations from site to site, 
reflecting a distribution of local environments "seen" by the electrons. 

Finally, in order to obtain a closed set of DMF self-consistency conditions, we 
need to specify a procedure that relates the lattice Green's functions Gj k to the 
solution of the local dynamical problem, as defined by the local effective action 
S e ff(i). We first note that the local action S e ff(i) is identical to that describing 
an Anderson impurity (AI) model (Anderson 1961, Hewson 1993) embedded in 
a sea of conduction electrons described by a hybridization function Aj(a> n ). The 
solution of this AI model then uniquely defines the corresponding self energy £j , 
which is given by 

Ei(w n ) = i0J n + (1 - Si - AiiuJn) - (G^K)) -1 , (5.7) 
where the local Green's function 

G^ C K) =< c\{u n )ci{u n ) > loc , (5.8) 
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is evaluated with respect to the local effective action S e ff(i). 

Next, we follow an "exact eigenstate" strategy, and define the "bare" lattice 
Green's functions G?- as the exact lattice Green's functions evaluated for the same 
realization of disorder {£»}, in absence of interactions. The "full", i. e. interaction- 
renormalized lattice Green's functions are within statistical DMF theory then 
defined by 

G ij = G° ij [£ i ^e i + X i (Lu n )}, (5.9) 

closing the set of DMF self-consistency conditions. 

We emphasize that a similar relationship relates the exact interaction-renormalized 
lattice Green's functions to their noninteracting counterparts. However, in the ex- 
act formulation, the self energies describing the interaction renormalizations are 
in general nonlocal in space, as well as in frequency. Having this in mind, we 
can describe the statistical DMF theory as a requirement for these interaction 
dependent self energies to assume a strictly local character. 

For an arbitrary lattice, an iterative procedure for solving the above set of 
statistical DMF equations could be obtained by 

(i) making an initial guess for the form of Aj(a> n ), 

(ii) solving the corresponding AI models on every lattice site, 

(iii) use the resulting Sj-s to calculate the full lattice Green's functions, 

(iv) calculate the new values of Aj(u; n ) and go back to step (ii). 

This procedure can be carried out for an arbitrary lattice in any dimension, 
but the procedure can be fairly time consuming due to the need to compute all 
the off-diagonal components of the lattice Green's functions Gij. However, similar 
computations have already been carried out to study the interplay of correlations 
and disorder in a Hartree-Fock (HF) approach (Yang & MacDonald 1993). In- 
terestingly, such a lattice HF treatment for a fixed disorder realization can be 
obtained as a further simplification of the statistical DMF theory, if the local 
Anderson impurity models are themselves solved at the same level. Of course, 
an exact solution of the statistical DMF equations goes well beyond such a HF 
treatment, since it can describe inelastic processes, which are expected to be 
of particular importance at finite temperature and in metallic non Fermi liquid 
phases. 

In practical terms, it is important to identify specific models where the solution 
of the statistical DMF equations can be simplified. The situation is particularly 
simple in the case of a Bethe lattice (Cayley tree). Because the absence of loops 

on such lattices, only local Green's functions G^{u n ) appear in the expression 

for Ai(uj n ). Furthermore, in this case the objects G^{u n ) can be computed from 
a local action of the form identical as in Eq. (5.2), except that in the expression 
for Aj(uj n ), the sum now runs over z — \ neighbors, excluding the site i. We thus 

conclude that the objects G^(o> n ) are related by a stochastic recursion relation, 
that involves solving Anderson impurity models with random on-site energies £j. 
In the non-interacting limit, the recursion relations can be written in close form, 
and reduce to the exact solution of the disordered electrons on a Bethe lattice 
(Abou-Chacra et al. 1971). 
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6. Order parameters for the metal to insulator transition: the 

disordered case 

In attempting to describe any phase transition, a crucial step is to identify 
appropriate order parameters that can characterize the qualitative differences 
between the different phases of the system. 

Since the basic focus of the dynamical mean field approach is the one particle 
Green's function, it is useful to write down an explicit expression for the Greens 
function of a lattice system with a one body Hamiltonian Hij and a general two 
body interaction term in terms of a self energy Sjj(icj). 

P _1 K)]»i = [tijiVn ~ Hij - T,ij(uj n )]. (6.1) 

Here, we have used a matrix notation, for a fixed realization of disorder, so that 
the self energy £^(0;) describes the interaction induced renormalizations of the 
Green's function. A quasiparticle picture emerges under a fairly general assump- 
tions of regularity of the interaction self energy at low frequencies. In this case 

EtfH « Eij(0) + ([Z-% - IK + 0{uo 2 n ). (6.2) 

and the Greens function at low frequencies assumes a very transparent form 

G(oj n ) = ^ZG QP {uj n )^Z (6.3) 

with the quasiparticle Green's function defined as 

[G^]- 1 =u n I- s/Z(H + X(0))VZ. (6.4) 

Introducing the low energy eigenvalues of the quasiparticle Hamiltonian E n and 
their corresponding eigenvector \n >, we can write 

, ^ < x\\fZ\n >< n\\[Z\x' > 
G(x, x ) « ^ j- — . (6.5) 

So far we have just repeated the derivation of the disordered Fermi liquid frame- 
work of dirty metals ( Castellani et al. 1987). In that reference those concepts 
were used to interpret the field theoretical coupling constants of Finkelshtein's 
non linear sigma model in terms of Fermi liquid parameters. Implicit in that 
framework was the assumption that the spatial fluctuations, or the sample to 
sample fluctuations were not too large. This assumption is valid for weak disor- 
der. The main advance of the dynamical mean field approach is the ability to 
deal with strong spatial fluctuations. 

The results of the next section indicate very different behavior for typical and 
average quantities, reflecting large spatial fluctuations. These differences may 
be the source of the finite scale divergences that were encountered in the field 
theoretical approach (Finkhelstein 1983, 1984) which, from the very begining, 
carries out an average over the disorder. 

We now use the previous developments to introduce the order parameters which 
are relevant to the DMF theory. We begin with defining the local density of states 
at zero frequency 

Pi = —ImGu(p+), (6.6) 

7T 
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and the corresponding local quasiparticle density of states 

p QP = _I /mG <f (0+ ). (6.7) 

Their ratio defines the second order parameter 

Zi = -gp- (6.8) 
Pi 

which could be dubbed the local quasiparticle weight. In the DMF framework, 
all the previous equations simplify because the self energy is local (Ejj(cj n ) = 
SijT,i(uj n )) and Zi can directly be calculated from the expression 

d i- 1 



Zi = 



1 - ^i2e[S i (w„)]| WTl =o+ 



(6.9) 



When these parameters are mostly uniform, i.e. site independent their averages 
are precisely two of the coupling constants in the field theoretical approach of 
Finkelshtein. In the statistical mean field approach we focus instead on whole 
distributions (or on typical values) of these parameters. 

We now turn to the physical interpretation of the order parameters, or more 
precisely their statistical distribution. The low energy physics is described by two 
parameters, pi and Zi, which are associated with the height and the width of 
the resonance in the spectra of the local impurity problem. Physically, p, L can be 
interpreted as a the density of states for adding or removing and electron from 
a specific site; Zi is related to the energy scale (or timescale) over which the 
quasiparticle picture applies. On the metallic side of the transition, the electron 
behaves as a local magnetic moment up to a time scale proportional Z~ x , but as 
a delocalized quasiparticle over longer timescales. The vanishing of Zi signals the 
conversion of quasiparticle degrees of freedom into local moments, which takes 
place at the Mott transition. 

To gain further insight into the physical content of these order parameters, we 
concentrate on transport properties. The fundamental difference between a metal 
and an insulator is defined by the ability of the electron to leave a given lattice 
site - to delocalize. In early work (Anderson 1958) Anderson pointed out that, at 
least for noninteracting electrons, this property could be examined by evaluating 
the lifetime of an electron on a given site. Very generally, the inverse lifetime 
simply measures the width of the the local resonance level. We thus expect 

- ~ Im(G ii (0+))- 1 . (6.10) 

Ti 

Within DMF, the local Green's function takes the form 

Gu(uj n ) = [iu n + n-Ei- Ai(uj n ) - £j(o; n )] _1 , (6.11) 

and since for Anderson impurity models (see for example Hewson 1993) Im Ej(0+) - 
0, we conclude that the desired local lifetime can be directly related to the local 
hybridization function as 

-~ImAj(0+). (6.12) 

n 

We thus expect ImAj(O) to vanish whenever the system is insulating, but to 
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remain finite in a metallic regime. We emphasizing that the same qualitative 
behavior should be expected for the local density of states (LDOS) pi as can be 
readily seen from Eq. (5.3). We can thus use the LDOS as an order parameter that 
discriminates a metal from an insulator. In a random system, pi will fluctuate 
from site to site, and we need a whole distribution function to fully characterize 
the approach to the transition. 

In particular, in the Mott insulator, there is a "hard" gap of order U on every 
lattice site, so even the average DOS discriminates the Mott insulator from the 
metal. The situation is more complex as the Anderson insulator is approached. 
Here, the local spectrum is composed of a few delta-functions (discrete, bound 
states), separated by gaps, but the average DOS remains finite. In contrast to 
the Mott insulator, in the Anderson insulator the sizes and positions of the local 
gaps fluctuate, but in both typical site has a gap at the Fermi energy. A 

natural order parameter is therefore the typical DOS, that is represented by the 
geometric average 

Ptyp =< P >geom= 6Xp{< In pi >}. (6.13) 

This quantity is found to vanish at the Anderson transition, in contrast to the 
average DOS, which is not critical. 

On the metallic side of the transition, the distribution function of a second 
quantity, the local quasiparticle (QP) weight, is necessary to characterized the 
low energy behavior near the transition. Important information is obtained from 
the typical value of the random variable Zi , defined as 

Z typ = exp{< In Zi >}, (6.14) 

which emerges as a natural order parameter from previous studies of the Mott 
transition. 

Finally, we define the averaged QP DOS by 

pgr =< p? p > ■ (6.i5) 

This object is very important for thermodynamics, since it is directly related 
to quantities such as the specific heat coefficient 7 = C/T, or the local spin 
susceptibility xioc- Note that in absence of interactions, p^ reduces to the usual 
(algebraic) average DOS, which is not critical at a U = Anderson transition, 
but it is strongly enhanced in the vicinity of the Mott transition. 

It is instructive to discuss the behavior of these order parameters in the previ- 
ously studied limiting cases. In the limit of large lattice coordination, the spatial 
fluctuations of the bath function Ai(uj n ) are unimportant, and there is no qual- 
itative difference between typical and average quantities. In the Mott insulating 
phase of a periodic solid, there is a gap in the density of states, while there is 
a finite density of states on the metallic side of the transition. As the MIT is 
approached from the metallic side, pt yp remains finite, but Zt yp is found ( Do- 
brosavljevic & Kotliar 1993,1994) to linearly go to zero. 

Another well studied limit is that of noninteracting electrons on the Bethe 
lattice, which is known (Abou-Chacra et al. 1971, Efetov 1987, Mirlin &; Fyodorov 
1991) to display an Anderson transition. The average DOS is finite both in the 
insulating and in the metallic phase, and is non critical at the transition. Similarly, 
by definition Zt yp = 1 in this noninteracting limit, so it also remains non critical. 
On the other hand, the typical density of states pt yp is finite in the metal and 
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zero in the Anderson insulator. This quantity is critical, and is found to vanishes 
exponentially with the distance to the transition. 

The definitions of the the order parameters that we have proposed are not 
restricted to the statistical DMF framework, which is simply used as a specific 
calculational scheme. The same definitions can in principle be used in other ap- 
proaches that can calculate local unaveraged values of the local DOS pi(ui n ) or the 
local part of the frequency dependent self energy T,a(uj n ) due to the interactions. 

7. Results 

The dynamical mean field theory maps the insulating phase of the model onto 
a collection of Anderson impurity models each one of them is embedded in an 
insulating bath. An Anderson impurity in an insulator, away from particle hole 
symmetry, can either have a doublet ground state when the coupling to the 
enviroment is weak or a singlet ground state when the coupling to the bath 
exceeds a critical value. In this strong coupling limit (Dobrosavljevic & Kotliar 
19926), a bound state is pulled from the continuum formed by the bands of the 
insulator. From these general consideration, we obtain a two fluid picture of 
the insulating phase: there are sites which have a local moment down to zero 
temperature while other sites quenched their spin by exchanging it with a strongly 
coupled neighbor. 

The number of sites in a doublet state, in the insulating phase of the system 
form "Mott droplets" . Based on our experience on the Mott transition in clean 
systems, we expect a larger density of states when the system undergoes a metal 
to insulator transition into an insulating phase with a large number of doublets. 
This result in at least two different regimes depending on whether the insulating 
phase has a high or low concentration of sites in doublet states. 

We have considered a z = 3 (Dobrosavljevic & Kotliar 1997). Bethe lattice, 
in the limit of infinite on-site repulsion U at T = and fixed average density n 
in the presence of a uniform distribution of random site energies £j of width W. 
To calculate the probability distributions of pj and Zj we used a simulation ap- 
proach, where the probability distribution for the stochastic quantity (co n ) is 
sampled from an ensemble of N sites, as originally suggested by Abou-Chacra et 
al. (1971). In order to solve Anderson impurity models for given bath functions 
Aj(uj n ) we use the slave boson (SB) mean-field theory (Barnes 1977, Read &i 
Newns 1983, Coleman 1987), which is known to be qualitatively and even semi- 
quantitatively correct at low temperature and at low energies. In agreement with 
heuristic arguments, we expect the results to be a strong function of the density 
n. In order to illustrate this behavior, we have carried out explicit calculations 
for both low electron density n = 0.3, and high electron density n = 0.7 (i.e. 
close to half filling). We emphasize that in the clean limit (W = 0) the behavior 
is qualitatively identical for the two values of the density, and the system remains 
metallic, with only the value of the effective quasiparticle mass m* ~ Z~ x be- 
ing a function of n. Qualitatively different behavior is found as the disorder is 
introduced. 

We first describe the evolution of the probability distribution of the local quasi- 
particle weights Zi, as the disorder is increased. The sites with Z\ <C 1 represent 
(Milovanovic et al. 1989, Dobrosavljevic & Kotliar 1993,1994) disorder-induced 
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Figure 1. DOS statistics in the low filling (high doping) regime. The full distributions P(ln p) are 
presented for increasing amounts of disorder. We find that the maximum, i. e. < lnp > shifts, 
as the transition is approached. Note also the extremely large width of the distribution, so that 
p now spans many orders of magnitude. 

local magnetic moments, and as such will dominate the thermodynamic response 
(see the definition of Pqp)- For weak disorder we expect relatively few local mo- 
ments and the quasiparticle weight distribution is peaked at a finite value. As 
the disorder is increased, the distribution of Zj-s broadens. At a critical value of 
the disorder W = W n fi, the form of this distribution assumes a singular form 
(Bhatt k Fisher 1992, Dobrosavljevic et al. 1992a, Miranda et al. 1996, 1997) , 
leading to anomalous thermodynamic response characterized by a diverging mag- 
netic susceptibility x an d specific heat coefficient 7. This behavior was found both 
for high and low density, in remarkable agreement with experiments carried out 
on for uncompensated and compensated doped semiconductors. Interestingly, a 
similar transition to a non Fermi liquid metallic phase, well before the MIT, has 
been found from the field-theoretical approaches in 2 + e dimensions (Finkelshtein 
1983,1984; for a review of subsequent work on this issue, see the review by Belitz 
k Kirkpatrick (1994)). 

As the level of disorder is increased further, a metal insulator transition is 
reached at the second critical value of disorder W = W c . At this point, the 
typical quasiparticle weight Zt yp vanishes, a behavior reminiscent of the clean 
Mott transition. While this behavior is found at both high and low density, the 
examination of the second order parameter pi reveals striking density dependence 
in the vicinity of the transition. The difference are most clearly displayed by 
plotting the evolution of the distribution function P(pi) as a function of doping. 
At low filling n = 0.3, we find that the with of this distribution becomes extremely 
large, spanning many decades, reflecting huge spatial fluctuations of the electronic 
wave- functions, as shown in Fig. 1. At the same time, the most probable value 
Ptyp is found to dramatically decrease, vanishing linearly with the distance to 
the transition. In contrast to this result, a very different behavior is found at 
high density n = 0.7. As we can see from Fig. 2, in this case the distribution 
width again broadens, albeit in a somewhat slower fashion, as the transition is 
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Figure 2. DOS statistics in the high filling (low doping) regime. W show the probability distri- 
bution for In p, as a function of disorder. We clearly see that while the distribution broadens, 
< \nq >« cons. 

approached. However, in this case the most probable value pt yp is only weakly 
modified with increasing disorder, and is found to approach a finite value at the 
transition. Recalling that the statistics of pi are intimately related to transport, 
our results strongly suggest that kinetic coefficients such as the conductivity will 
strongly depend electron filling, as the MIT is approached. These expectations 
should be confirmed by explicitly calculating the critical behavior conductivity 
at different values of the density. 

Finally, we mention an important feature of the insulating state in presence 
of both the strong electronic correlation and disorder. By explicit calculations of 
both the average and the typical density of states we have demonstrated that 
the introduction of disorder, fills the Mott gap with localized electronic states. 
The average DOS is therefore finite in the insulator. Nevertheless, the MIT still 
retains a definite Mott character, with a finite fraction of electrons turning into 
localized magnetic moments, and the typical quasiparticle weight Z typ vanishing 
at the transition. Our results explicitly show that the metal insulator transition in 
presence of both interactions and disorder is a qualitatively new type of transition, 
having well defined signatures of both the Anderson and the Mott route to 
localization. 



The metal to insulator transition continues to be one of the central problem 
in condensed matter physics. Building on the fundamental concepts introduced 
by Mott and Anderson there have been several attempts at the construction of 
a microscopic theory. We have summarized in this contribution some aspects of 
the dynamical mean field approach to this problem. 

Deeply rooted in local physics, it combines the chemical and physical aspects 
of the problem. We regard the construction of a mean field theory as an essential 
first step before a more complete treatment including Gaussian and non linear 



8. Conclusions 
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fluctuations is carried out. In the clan limit the DMF description of the Mott 
transition has already given insights into puzzling aspects of transition metal 
oxide physics. 

Since the approach emphasizes the local environment, it has a conceptually 
very simple extension to the disorder case: the statistical dynamical mean field 
theory. The theory explicitly incorporates both the Mott and the Anderson route 
to localization, and thus provides a consistent description of the transition, in- 
terpolating between the respective limits of no disorder and no interaction. This 
key feature seems to be missing in a recent field theoretical formulation of the 
Anderson- Mott transition (Kirkpatrick & Belitz 1994, 1995). 

The statistical mean field theory reproduces many remarkable features of doped 
semiconductors. A transition to a Griffiths phase (Bhatt and Fisher 1992, Do- 
brosavlevic et. al. 1992, Lakner et. al. 1994) where local moments coexist with 
conduction electrons down to arbitrary low temperature precedes the true metal 
to insulator transition. The local properties such as the density of states depend 
in a distinct fashion on the level of doping (compensation). The metal to insulator 
transition in the low doping region (uncompensated) has a larger typical density 
of states and is closer to a Mott transition than the corresponding transition 
at larger doping levels (compensated situation). The latter contains ingredients 
from both the Anderson and the Mott transition. 

There are several issues that deserve further investigation. A detailed study 
of the transport properties has to be carried out. One should also elucidate the 
interplay of the short wavelength fluctuations with long wavelength modes. The 
latter are presumably described by the non linear sigma model approaches that 
have been investigated near two and six dimensions. A second element still missing 
form the mean field theory is the inclusion of short range magnetic correlations. 
In our view these are the most pressing problems that need to be address on the 
road to a comprehensive theory of the Mott- Anderson transition. 
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